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The probability distribution of the reflection coefhcient for light reflected from a one-dimensional 
random amplifying medium with cross- correlated spatial disorder in the real and the imaginary parts 
{^JQ. of the refractive index is derived using the method of invariant imbedding. The statistics of fluc- 

^ ' tuations have been obtained for both the correlated telegraph noise and the Gaussian white-noise 

I models for the disorder. In both cases, an enhanced backscattering (super-reflection with reflection 

coefhcient greater than unity) results because of coherent feedback due to Anderson localization 
and coherent amplification in the medium. The results show that the effect of randomness in the 
imaginary part of the refractive index on localization and super-refiection is qualitatively different. 



^ ' Light propagation in disordexed media and the associated Anderson localization of a wave in both the active as well 
Ph . as the passive random mediaEmhas been studied extensively. In recent years, there has been increased interest in light 
' propagation and lasing in active random media supported by the several experiments carried out on these systemstl^Ll. 
However the experimental findings of a narrowed spectral emissionO'EI and a pulse narrowing of the emissionQ'u above 
a well defined threshold of pumping could be explained merely as an effect of the loHg-jdiffusive pathlengths in a ran- 
, dom medium with gain and the consequent amplified spontaneous emission (ASEpM. Mare recently, the observed 

■ super-narrowing of the emitted spectra fcam strongly scattering semiconducting powdeiO and from weak scatter- 
I ] ers dispersed in high gain organic mediat3 has been attributed to coherent feedback caused by recurrent multiple 

'"O ' scatteringEj. It is, however, still debatable if the wave amplification is due to the predicted synergy between wave 
^ confinement by Anderson localization and coherent amplificationa. 

^ ■ In all these studies, the active random medium is considered to scatter the propagating wave (light) due to fluc- 
tuations in the real part of the refractive index {rjr) (real potentials) while the coherent amplification is modelled by 
a phenomenological spatially constant imaginary part of the refractive index (77^). However, it would be of interest 
' to examine the effect of a spatially fluctuating imaginary part of the refractive index as well. More so, as the scat- 
j tering micro-particles (e.g. polystyrene microspheres, Ti02 rutile particles) used in the experiments are not active, a 
corresponding mismatch in the imaginary part of the refractive index is found to exist. It has been pointed out by 
QQ Rubio and KumaiHj that a mismatch in the imaginary part of the refractive index (imaginary potential) would always 
cause a concomitant reflection (scattering) in addition to the absorption or amplification. Mismatch in r]i alone in an 
0^ , amplifying medium (negative imaginary potential) with no mismatch in rjr can cause resonant enhancement of the 

■ scattering coefficients. In fact, the reflection and the transmission coefficients can even diverge as can be seen from 
' the simple example of a single imaginary 6 potential in one dimension. This would correspond to the experimen- 

tal situation where the scatterers (polystyrene microspheres, say) are suspended in a fluid with the same rjr (index 
matching fluid) in which a laser dye is dissolved and optically pumped. The scattering caused by the fluctuations in 
, rji would, therefore be expected to have non-trivial effects on the wave propagation in the medium. 

Q , The transmittance across a randomly amplifying and absorbing chain was recently considered by Serlli numerically 
O ' and was shown to decay exponentially with the increase in length of the chain, presumably due to localization. But 
the effects of the fluctuation in the imaginary part of the refractive index on lasing in such random media has not 
been studied so far. In this work, we consider the statistics of the non-self- averaging fluctuations of the reflection 
coefficient for light incident on a one-dimensional active random medium with spatial correlated disorder in both the 
imaginary part as well as the real part of the refractive index. 



- ' — \ 

X 



We consider a one-dimensional active disordered medium of length L with a random complex refractive index 
r],0 < X < L. For simplicity, polarization effects are neglected and light is assumed to be a scalar wave. A physical 
realization of interest here would be an Er^~^ doped and pumped, polarization maintaining optical fibre intentionally 
disordered along its length. Further, only the linear case of the gain/ absorption being independent of the wave 
amplitude is considered and non-linear features such as gain saturation are not considered. Here we would like 
to emphasize, that our treatment is for the possibility of super-reflection (r > 1) i.e., for an amplifier and not an 
oscillatorE3EII. The complex wave amplitude E(x) obeys the Helmhotz equation inside the medium. 
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where k is the wave vector m the medium (/c^ = w^/c^eo) and ri{x) — rjr{x) + i[f]i + r^i^x)] is the complex refractive 
index. Here rjr{x) and rji{x) are random and fji is a constant representing the average amphfication or absorption in the 
medium according as f]i is negative or positive. It is well known that Eq.(|^) can be transformed to give an equation for 
the evolution of the emergent quantity, namely, the complex amplitude reflectipiixoefficient R{L) — [r(L)]^/^ exp[iO{L)] 
as a function of the sample length L, via the method of invariant imbeddingESllj as 



dR{L) 
dL 
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Equation. (||) is a stochastic differential equation and we are interested in the corresponding Fokker-Planck equation 
for the probability distribution P{r, 9; L) which can be readily obtained following the standard procedures. Thus, 
let n(r, <kJ^) be the density of points in the. (r, 0) phase space. Now n(r, 6; L) must satisfy the Stochastic Liouville 
equatiortj, and by the van Kampen lemmaE3, the probability distribution function P{r, 9; L) = (n(r, 9\ L))n,..rii , where 
the angular brackets denote averaging over all the realizations of the random refractive indices rjr and ry^. 

A. The Gaussian S correlated (white-noise) disorder 

First, let us consider the simplest case namely that of a Gaussian ^-correlated (white-noise) model. In this model, 
Tjr and r]i are assumed to have S correlated Gaussian distributions with {rir{L)) — 0, {r]i{L)) — 0, {•qr{L)rjr{L')) = 
A'^S{L — L') and {rii{L)r]i{L')) — AfS{L — L'). This model would most appropriately describe the case of a continuous 
random medium such as a laser-dye dopied gel or intralipid suspensiono, where the fluctuations in rjr and rji are 
uncorrelated. Using the Novikov theoremE3 to average over all configurations of rjr and ry^, we obtain in the Random 
Phase Approximation (RPA) {i.e., P{r,9) = P(r)/27r), 
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where the linear operators Lr. and Lj are given by 
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and the non-dimensional sample length I = 
Af/max{Al,Af} and A = 27]i/max{Aj, Af}k 
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(j)r = Al/max{Al,Aj}, (t>r = 
is the amplification length in the 



medium defined by the average of the imaginary part of the refractive index and max implies the superior value of 
the arguments. The RPA is known to be valid in the the weak disorder limit, klc >> 1, where Ic is the localization 
lengthEfl. We point out that even if rjr and rji were cross-correlated, the final equations do not differ in the RPA for 
the white-noise model (because {'Li'L2P)e = see equations (9), (10)). 

The asymptotic I —>■ 00 limiting solution of Equation. (|^) obtained by setting dP/ dl = is given by, 



P(r; 00) = Po 
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+ (t)i)r + b4)i - (t>r - ^ 



+ (j)i)r + 50i -0^ + 7 



-A/7 



> 20, 
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where 7 = ■y/l20i|0r — 20^1 and Pg is a normalization constant given by [j^ P{r, oo)dr] ^ . The limit I 00 implies 
physically L >> Ic. This expression goes over straightforwardly to the result of Pradhan and Kumaro in the limiting 
case of pure real disorder (0; = 0). Thus the statistics qualitatively differ in the two regimes for an amplifying 
medium : (i) when the real part disorder dominates (0^ > 20^) and (ii) when the Imaginary part disorder dominates 
(0, < 20/). 
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We have also solved equation(|^ numerically for finite length to investigate the approach to the asymptotic forms 
given by eq.(||). In Fig.(l), the plots of P(r, I) for the case of real disorder dominating {(pr > 2(pi) for different lengths 
of the medium are shown. The probability distribution for the case of a pure imaginary mismatch {(pr = 0), with the 
real part 77^ being index- matched is shown in Fig. (2). The line joining the dots in both the figures corresponds to the 
asymptotic P(r; 00) solution. In the case of amplifying medium, the value of reflectivity (rmax) at which P(r; /) peaks 
increases with the average value of the amplification factor |^|. For the case of imaginary part disorder dominating, 
P{r; I) has a peak at small values of the reflectivity even for moderate values of the amplification. In the case of an 
absorbing medium with the imaginary diorder dominating, the probability distribution has a monotonic decreasing 
behaviour and is maximum at r = 0. A finite probability of reflection at r > 1 in the absorbing case and at r < 1 
in the amplifying case {A < 0) is recognised to be a consequence of the two-sidedness of the white-noise process for 
the complex refractive index, which allows the imaginary part of the refractive index {ffi + r]i) to take on locally both 
positive and negative values for any given value of the average. It should be noted that this limiting form of P{r, 00) 
gives a weak logarithmic divergence for (r) (for ^ ), regardless of the sign of A for both absorption and amplifica- 
tion. Thus amplification has a much more drastic effect on the reflectivity than attenuation. The white-noise process 
allows the local fluctuations in rji to be very large and the effect of a finite mean value fji is small. It is thus a case of 
the fluctuations dominating over the mean. We also find that the numerical solutions saturate to the limiting forms 
for 1^1. So most of the reflection occurs fpfjtf* within a localization length. This enhanced backscattering is quite 
different from that caused by light diffusionErEJ. In the latter case, the distribution of optical path length, because of 
exponential growth of wave amplitude due to coherent amplification in one-dimension, gives P^(r; 00) ~ ln(r)^/^/r 
for r >> 1. This decays much slower than the P{r; 00) for r — > cxd, as given by eq.(p|). 



B. Correlated telegraph disorder 



In the case of the white-noise disorder, the imaginary part of the refractive index was allowed to take on both 
positive and negative values i.e., the medium could be locally both amplifying or absorbing. With a view to studying 
purely amplifying/absorbing random media, we use the telegraph disorder model to describe the fluctutations in the 
refractive index. Moreover, since the gain/absorption coefficient is physically always bounded from above, the fluc- 
tuations in the imaginary part of the refractive index are better described by this dichotomic Markov process {i.e., 
spatial Telegraph noise). Further, we recognize that in discrete random media such as microparticles suspended in 
a laser dye solution used in experiments, the real and the imaginary parts of the refractive index fluctuate spatially 
in the same manner and can, therefore, be described by the same stochastic process. A telegraph noise with a finite 
correlation length is most appropriate to describe such a situation. Accordingly, we will choose rir{L) = ax{L) and 
rji{L) = Px{L) with an average value for the imaginary part fji. Here x(^) is taken to be a dichotomic markov process 
which can take on the values ±x such that {x{L)) = and {x{L)x{L')) = exp (— Fji/ — where F^^ is the 
correlation length in the medium. 

Now, defining as before, P{r,9;L) ==jMl{r,6; L))^ and W{r,9;L) = {x{L)Il{r,9; L))^, and using the "formulae of 
differentiation" of Shapiro and LoginovEJ to average over the dichotomous configurations of x(i), we obtain 
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where the linear operators Li and L2 are 
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1.2= k 
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We thus get a closed system of equations for P{r,9,L) and W{r,9, L). These equations go over correctly to the 
corresponding eq.(|^) in the white- noise limit obtained by taking the limit x^ — * 00, T — 00 while keeping x^/r = 
constant. In this limit, the equation for P{r,d;L) becomes autonomous i.e., it gets decoupled from W{r,9; L). 



In the RPA and in the asymptotic limit L —t 00, these equations simplify to 
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where Lr and Lj are given by eq.(^ and eq.(|^) and A = 2Tfji/x^. Interestingly in the case of the pure real part 
disorder (/3 = 0), the form of the telegraph noise equation for P(r; oo) is identical to that for the white- noise case, 
but with the coefficient A = 2Tr]i/kx'^- Similiarly, in the case of the pure imaginary part disorder (a = 0), the form 
of the telegraph noise equation for P{r; oo) is again identical to that for the white-noise case, but with the coefficient 
A = 2Tfji/k[x^ — fji). However, for Px < IVil^ the imaginary part of the refractive index is always positive (absorbing) 
or negative (amplifying). Hence the solution for these two cases is also given by eq.(H), the solutions being valid in 
the interval < r < 1 for the absorbing medium, and 1 < r < oo for the amplifying medium. Outside the intervals, 
the probability density P{r;L) vanishes. 



A complete solution for the eq.(|ll|)and eq.(|12[) is obtained as 

1 1 



P(r; oo) = Po 
J±(r) = - 



e+(i + c+r + r2) e-(i + c-?- + ^^) 



exp[-2^(/+(r)+/_(r))] 



(13) 



1 



In 



.(2) 



: tan 



.(1) 



C ± + 2r 



IC±I >2 



IC±I <2 



where ^± = a 
.1/2[C± 



^l/2 



(1) 



-1/2[C± 



4) 



1/21 



— 4)^' 2] and Pq is a normalization coefficient. This expressions become the same as given 



(2) 

by eqn.(||) in the white noise limit (x^ ^ 0, F ^ oo and x^/r being constant). 

The solutions for one-sided disorder in the imaginary part exhibit three qualitatively different behaviours corre- 
sponding to choices of the parameters a, /3 and fJi (x is an arbitrary constant and can be set to unity without loss 
of generality). First, we note that the case of + P'^ — (3f]i/x = 0; corresponds to a singular perturbation of the 
differential equation for P(r; oo). This condition can be interpreted as a threshold condition by noting that the local- 
ization length is given hy ~ (a^ -1-/3^) and the effective amplification length is given by l^^l^p — (3-qi. This condition 
then corresponds to a matching of length scales in the problem, = lamp- In the regime where the amplification 
dominates the localization (a^ -t- 0^ 
in the region of interest (1 < r < oo) 
appreciably, as can be seen from Fig. (3a). For [a^ + 0'^ — Pfji > or Ic < lamp), a natural boundary arises for the 

('2) 

solutions of the equation at r_ which falls in the domain of physical interest (1 < r < oo). Now the solutions given 
by the expression (p^) are valid in the range r'^'' < r < oo with P{r; oo) = outside. In this regime the localization 
dominates {la < lamp), if 2^/[e-(C- - ^)-^'^] 



' Prji < or Ic > lamp), the solutions exhibit a monotonic decreasing behaviour 
Here the disorder in the real part (a) is small and does not affect the statistics 



(2) 



P(rL^^ oo) 



(Fig. (3b)). The value of ' 



> 1 and we have a broad distribution with peak at Vmax > t 
is large for small disorder in the real part (a^ + P'^ — p-fji > 0) 



and 
and 



decreases as a increases. The behaviour in this region is dominated by the disorder in the real part of the refractive 
index. A third qualitatively different behaviour occurs for^c < lamp and 2 A/ [^_(C^- 4) -1/2] < 1. Then the expression 

given by eqn.(|l3|) diverges at r^_^\ This divergence is, however, normalizable implying that P(r; oo) is peaked (in fact, 
sharply) at that point. This behaviour can be readily understood by noting that the second condition which can be 
rewritten as ■qi^iV /k)'^ < 3/3(|?7i| — ^la^ + 2P{\'fii\ —/?)], is basically a condition on the correlation length [Icorr = T~^). 
This condition is satisfied for small F (large Icorr)- Then the reflection is essentially from a single potential barrier 
and thus has a sharply defined value. It should be noted that, as a ^ oo, P(r; oo) 6{r — 1), as expected. 

The solutions for the case of a two-sided disorder for the imaginary part (/3 > |?7i|) are similiar to the solutions 

(2) 

for the white noise case. It should be noted that there does not exist real r_ ' which falls into the physical region of 
interest (0 < r < oo). In this case the large disorder in the imaginary part (/3) causes the effects of localization to 
dominate. However, in all cases of amplification, for a finite A and + P'^ — Pfji ^ 0, there is a universal l/r^ tail for 
the P(r; oo). For the case of pure imaginary disorder(a = 0), we similiarly see a monotonically decreasing behaviour 
of P(r; oo) with r for one-sided disorder {P < \f]i\ or Ic > lamp) (Fig-(3c)), and a P(r; oo) with a peak for two-sided 
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disorder (/3 > |77i| or Ic < lamp) (Fig. (3d)). With increase in f3 for two-sided disorder, the epak shifts to smaller 
values of reflectivity as the effects of absorption show up, until for large enough /3, the peaks occurs at r = and we 
again have a monotonically decreasing P(r; oo). It should be mentioned that all these effects are seen for the case of 
absorption also, with the roles of rL^"* and r'^^' interchanged. 

In conclusion, we have studied the statistics of super-reflection from a one-dimensional disordered system with 
spatial randomness both in the real and the imaginary parts of the complex refractive index. We have discussed the 
models of disorder qualitatively applicable to experimental systems such as intentionally disordered optical fibres with 
gain (i?r^+ -doped) and obtained the probability distribution function of the reflectivity for the cases of a white-noise 
disorder and a correlated telegraph disorder. In both cases, an enhanced reflection results because of coherent feed- 
back due to Anderson localization and coherent amplification. In the case of white-noise disorder, the statistics are 
qualitatively different in the two regimes of the real part disorder dominating (A^ > 2Af ) and the imaginary part 
disorder dominating (A^ < 2A|). In the case of telegraph disorder, we obtain three qualitatively different behaviours 
for P{r; oo) depending on threshold conditions involving the localization length, the amplification length and the 
correlation length. Thus the fluctuation in the imaginary part of the refractive index is seen to have a non-trivial 
and qualitatively different effect on localization and lasing from such random media. Finally, as the phenomenon 
considered here is concerned with the issue of statistical fluctuations (noise) in a random amplifying medium, we 
propose for it the acronym RAMAN (Random Amplifying Medium And Noise). 
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FIGURE CAPTIONS 

Fig. 1. The probability distribution of reflectivity P(r; /) in the case of the white noise disorder given by eqn.(^ and 
the real disorder dominating (0^ = 1-0, (t>i = 0.1), for the different sample lengths indicated. The line joining the dots 
is the analtic result for P{r; oo). The amplification parameter is A = —0.25 . 

Fig. 2. The probability distribution P{r;l) in the case of the white noise disorder given by eqn.(|^) and a pure 
imaginary mismatch {(pr — 0) for different lengths of the sample. The line joining the dots is the analtic result for 
P(r; oo). The amplification parameter is A = — 1. 

Fig. 3. The probability distribution P(r; /) in the case of the correlated telegraph noise, (a) > lamp and (b) 
Ic < lamp are for one-sided disorder (/3 < \f]i\) with disorder in both the real and the imaginary parts, (c) Ic > lamp 
and (d) Ic < lamp are for two-sided disorder (/3 > \f]i\) and pure imaginary mismatch (a = 0). 
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(a) P = 0.3;ti = 1.0;r= 1.0 



(b) P = 0.5;ti = 1.0;r = 3.0 
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(c) a = 0.0;fi = 1.0;r= 1.0 
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